All-sky analysis of the general relativistic galaxy power spectrum 



o 

(N 
00 



Jaiyul Yoo 1: Q and Vincent Desjacques-fj] 

1 Institute for Theoretical Physics, University of Zurich, CH -8057 Zurich, Switzerland 
2 Lawrence Berkeley National Laboratory, University of California, Berkeley, CA 94720, U.S.A. and 
' ' Departement de Physique Theorique and Center for Astroparticle Physics, Universite de Geneve, CH-1211 Geneve, Switzerland 

(Dated: January 22, 2013) 

We perform an all-sky analysis of the general relativistic galaxy power spectrum using the well-developed 
spherical Fourier decomposition. Spherical Fourier analysis expresses the observed galaxy fluctuation in terms 
of the spherical harmonics and spherical Bessel functions that are angular and radial eigenfunctions of the 
Helmholtz equation, providing a natural orthogonal basis for all-sky analysis of the large-scale mode measure- 
ments. Accounting for all the relativistic effects in galaxy clustering, we compute the spherical power spectrum 
and its covariance matrix and compare it to the standard three-dimensional power spectrum to establish a con- 
nection. The spherical power spectrum recovers the three-dimensional power spectrum at each wavenumber k 
with its angular dependence fik encoded in angular multipole I, and the contributions of the line-of-sight projec- 
tion to galaxy clustering such as the gravitational lensing effect can be readily accommodated in the spherical 
Fourier analysis. A complete list of formulas for computing the relativistic spherical galaxy power spectrum is 
also presented. 

PACS numbers: 98.80.-k,98.65.-r,98.80Jk,98.62.Py 
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I. INTRODUCTION 

The past few decades have seen a rapid progress in large- 
scale galaxy surveys. The Sloan Digital Sky Survey (SDSS; 

and the Two degree Field Galaxy Redshift Survey (2dF- 
GRS; [2]) opened a new horizon in modern cosmology, map- 
ping three-dimensional positions of millions of galaxies. The 
Baryonic Oscillation Spectroscopic Survey (BOSS; H) and 
the Dark Energy Survey (DES; [4]) represent the current state- 
of-the-art galaxy surveys, together with the recently com- 
pleted WiggleZ Dark Energy Survey |5[. The exploration of 
the large scale structure of the Universe will continue with 
future galaxy surveys such as Euclid, 1 the BigBOSS, 2 the 
Large Synoptic Survey Telescope, 3 and the Wide-Field In- 
fraRed Survey Telescope. 4 As they will cover a substantial 
fraction of the entire sky and a wide range of redshift, these 
future surveys will measure galaxy clustering with stupendous 
statistical power, demanding thereby that the current theoreti- 
cal predictions be refined to achieve higher levels of accuracy. 

The galaxy distribution measured in galaxy surveys repre- 
sents a biased version of the underlying matter distribution. A 
traditional approach to analyzing galaxy number density fields 
is to utilize the power spectrum of its Fourier components. 
Since each Fourier mode evolves independently in the linear 
regime, the power spectrum measurements can be used to in- 
fer the primordial matter power spectrum and to extract cos- 
mological information (see for recent power spectrum 

measurements). While the power spectrum analysis merits its 
intuitively simple interpretation of the measurements in con- 
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nection with the underlying matter distribution, it assumes that 
the density fields are defined in a cubic volume. As the recent 
and forthcoming galaxy surveys cover a progressively larger 
fraction of the sky, the validity of the flat-sky approximation 
and the power spectrum analysis becomes questionable. 

At the same time, there exists a demand for large-scale mea- 
surements from the theoretical side. In the past few years the 
relativistic description of galaxy clustering has been devel- 
oped |9, 10]. The advance in theoretical development results 
from the finding that the observed quantities in galaxy cluster- 
ing such as the observed redshift and the galaxy position on 
the sky are different from quantities used to construct the ob- 
served galaxy fluctuation such as the background redshift and 
the unlensed galaxy position 111 111 . Those theoretical quantities 
are gauge-dependent, and the subtle difference in those quan- 
tities become substantial on large scales, where the relativistic 
effect becomes important |9lll0(]. The full relativistic formula 
of galaxy clustering can be analytically derived at the linear 
order in perturbations, providing a complete picture of galaxy 
clustering on large scales |9l [lOl Il2l4l6ll . Furthermore, it is 
shown [17] that these relativistic effects in galaxy clustering 
can be measured in future galaxy surveys, providing a great 
opportunity to test general relativity on cosmological scales. 

The relativistic formula for the observed galaxy fluctuation 
S° hs (z, x) is well-defined in observations, where it is a func- 
tion of the observed redshift z and angle x on the sky. Since 
the wide angle formalism Jl8l - l2lll has been developed to com- 
pute the two-point correlation function of the full Kaiser for- 
mula without adopting the distant-observer approximation, its 
extension to the full relativistic formula can be readily made 
|22[ l23ll . However, the resulting equation for the two-point 
correlation function is highly complicated, even for the full 
Kaiser formula, mainly due to the geometrical effect. Further- 
more, its physical interpretation is not as straightforward as in 
the power spectrum analysis owing to its nontrivial relation to 
the primordial matter power spectrum. 

Adopting the flat-sky approximation, Yoo [ 10] computed 
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the galaxy power spectrum, accounting for the relativistic ef- 
fects (see, also, ll6H ). While it is always possible to embed 
the observed sphere in a cubic volume with rectangular co- 
ordinates and to perform a power spectrum analysis, it be- 
comes difficult in principle to connect these large-scale mea- 
surements to the underlying theory, because the fiat-sky ap- 
proximation has a limited range of validity. However, it is 
shown [24] that, on these large scales where measurement 
uncertainties are significant, the systematic errors associated 
with the fiat-sky approximation are indeed negligible in the 
power spectrum analysis, if performed properly. This is con- 
sistent with the previous finding of liiclll7ll obtained by a sim- 
pler treatment. 

Here we present an alternative to the traditional power spec- 
trum analysis, based on radial and angular eigenfunctions of 
the Helmholtz equation. The spherical Fourier analysis has 
been well developed 1 25-27] in galaxy clustering, while its 
application was limited to the Kaiser formula. The observed 
galaxy fluctuation is decomposed in terms of Fourier modes 
and spherical harmonics, both of which provide a natural or- 
thogonal basis for all-sky analysis. In observational side, 
Fisher et al. [28] applied the spherical Fourier analysis to the 
IRAS Redshift Survey, and the method was further developed 
in Fisher et al. ll26tl to reconstruct the velocity and the gravita- 
tional potential fields. Theoretical refinement has been made 
in the past few years. Rassat and Refregier ll29ll computed 
the spherical power spectrum of the matter density with fo- 
cus on the baryonic acoustic oscillation feature, ignoring the 
redshift-space distortion (see OOll for the effect of the redshift- 
space distortion and nonlinearity). Shapiro et al. Bill used the 
spherical Fourier analysis of the redshift-space distortion and 
its cross-correlation with the CMB temperature anisotropics 
to constrain the growth rate of structure. The same technique 
is also applied to the weak gravitational lensing formalism 
in Heavens lf32ll to take advantage of the information on dis- 
tances to background source galaxies. Compared to the tradi- 
tional weak lensing, this spherical Fourier analysis is known 
as the 3D weak lensing, and it is shown in Kitching et al. f33h 
that the 2D tomography in weak lensing is just the 3D weak 
lensing with the Limber approximation. 

Drawing on this previous work, we perform a spherical 
Fourier analysis of the observed galaxy clustering. Account- 
ing for all the general relativistic effects in galaxy cluster- 
ing, we compute the spherical power spectrum and its co- 
variance matrix and compare them to the corresponding three- 
dimensional power spectrum. The organization of this paper is 
as follows. In Sec. [Til we briefly review the spherical Fourier 
analysis. We first discuss the basic formalism of the spheri- 
cal Fourier analysis in Sec. Ill Al and its application to galaxy 
redshift survey in Sec. MB I In Sec. [Till we present the full 
general relativistic description of galaxy clustering. We first 
describe the key equations behind the relativistic description 
in Sec. IHI Al and present its relation to the spherical Fourier 
analysis in Sec. IHI Bl In Sec. IHI Cl we provide the Limber ap- 
proximation to the spherical Fourier analysis of the relativistic 
description. The main results of the spherical Fourier analy- 
sis are presented in Sec. [IV] in which we discuss the spher- 
ical power spectrum of the relativistic formula in Sec. IIV Al 



and the measurement uncertainties associated with them in 
Sec. IIVBI Finally, we discuss the implications of our results 
in Sec. [V] Our numerical calculations are performed by as- 
suming a flat ACDM universe with cosmological parameters 
consistent with the WMAP7 results [34]. 



II. FORMALISM 

Spherical Fourier analysis of three-dimensional density 
fields has been well developed j25l - [27ll (see, also, ll35ll for a 
different approach). We briefly review the basics of spherical 
Fourier decomposition in Sec. Ill Al and discuss its application 
to galaxy surveys in Sec. Ill Bl accounting for issues associated 
with redshift distribution. 



A. Spherical Fourier Decomposition 

We begin by defining our normalization convention for 
Fourier decomposition. The position vector |x) and the 
Fourier mode |k) are normalized to satisfy the completeness 
relation: 

1 = J d 3 x|x)(x| = J -|^|k>(k|. (1) 

With the plane wave represented by 

(x|k) = exp(ik • x) , (2) 

the configuration and the Fourier space vectors are normalized 

as 

(x|x') =<J £) (x-x / ) , (k|k') = (27r) 3 ^(k-k') . (3) 

Given a scalar field <5(x) = (x.\6), its three-dimensional 
Fourier components are represented as 

5(k) = J d 3 xe- !k ' x 6{x) = J d 3 x (k|x)(x|<5) = (k|«5) , 

(4) 

and its ensemble average defines the power spectrum as 

(6(k)6* (k')> = P(k, k') = (2tt) 3 <5 d (k - k')P(fc) , (5) 

where the last equality holds if the power spectrum is rotation- 
ally and translationally invariant. 

With the normalization convention, we consider a complete 
radial and angular basis \klm) in a spherical Fourier space 
to decompose a scalar field in three-dimensional space. We 
define its representation in configuration space as 

(x.\klm) = J^kji(kr)Y lm (k), (6) 

where r — |x|, x = (6, <fi) is a unit directional vector of x, 
ji (kr) is a spherical Bessel function, and Y/ m (x) is a spherical 
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harmonics. The normalization coefficient is chosen such that 
the basis \klm) is orthonormal 

(k'l'm'\klm) = J d 3 x (k'l'm'\x){x\klm) 

= 5 D (k~k')6 U 'S mm , , (7) 

where we used 



drr ji(ar) ji(br) = — S (a ~ b) 
2ab 



(8) 



By expanding the plane wave in Eq. (0 

(x|k) = A-K^jtikrjY^&YUZ) 

lm 

= (27r) 3 / 2 ^^r^(k)(x|^m), (9) 

lm 

and using the completeness condition of the \klm) -basis 

(x|k) = J dk'^2{x\k'lm){k'lm\k) > (10) 

lm 

we derive the relation between our spherical Fourier basis and 
the usual Fourier mode 

(k\k'lm) = (2^) 3 / 2 t^ m (k) S D (k - k') . (11) 

Naturally, the spherical Fourier basis \klm) encodes the am- 
plitude k = |k| of the three-dimensional Fourier mode |k) and 
its angular direction Y/ m (k). 

Based on the |Mm)-basis, any scalar field in configuration 
space can be spherically decomposed as 

/>oo 

8(x) = (x\5) = / dk^2(x\klm)(klm\S) (12) 

J ° lm 

= J dk^2 k jl ( kr ) Y lm (x) Sim (k) , 



lm 



and its spherical Fourier mode <5/ m (fc) is related to the three- 
dimensional Fourier component as 



Sim(k) = (klm\S) 



i'k 



(2tt) 3 / 2 



d'kY t * m (k) S(k) . (13) 



Due to our normalization convention of the spherical Fourier 
basis, the spherical Fourier mode Si m (k) is slightly differ- 
ent from the usual coefficient of the angular decomposition, 
which is just the angular integral over k. Finally, the spherical 
power spectrum Si (k,k') is defined as 



(*, m (A0 s; m ,(k')) = 5 w s mm ,Si(k,k') 
i l {-i) l 'kk' 



(14) 



(2tt) 3 



d 2 kd 2 k' Y lm (k)Y l/m ,(k')(S(k)6*(k')) . 



For a rotationally and translationally invariant power spectrum 
in Eq. q), the spherical power spectrum is 

Si(k,k') =5 D (k-k')Si(k) =S D (k-k')P(k) , (15) 



and it reduces to the three-dimensional power spectrum 
Si(k) = P{k). Using the representation S(x) in configuration 
space, a similar but more convenient formula can be derived 
for the spherical Fourier mode as 



(16) 



S lm (k) - / d 3 x xl-kMkr) Y lm (x) S(x) 



and its spherical power spectrum is then 

2kk' 



Si(k,k') 



J d 3 Xl J d 3 x 2 ^( Xl )y; m (x 2 ) 

xji(kn)ji(k'r 2 ) (<J(xi)*(x 2 )) . (17) 



Fourier components and its power spectra are dimensionful: 



[5im(k)} - L 2 , 
[S(k)] = {P(k)] = {Si(k)}=L\ 
[Sti^k')] = L 4 , [P(k,k')] =L 6 . 



(18) 



B. Redshift Distribution and Survey Window Function 

In observation, we can only measure galaxies in the past 
light cone. Therefore, radial coordinates parametrized by the 
observed redshift carry special information, namely, time — 
quantities at two different radial coordinates are at two differ- 
ent redshifts. Furthermore, the mean number density n g (z) 
of galaxies evolves in time, such that their fluctuation field 
5 g should be properly weighted by their redshift distribution 
dfig/dz. This complication can be readily handled by decom- 
posing the number density field n 9 (x), instead of its fluc- 
tuation field 8 g {x) — n g (x)/n g (z) — 1, where the three- 
dimensional position vector x is understood as a function of 
its radial position r(z) — |x| and angular position x. 

Since the mean n g (z) is independent of angular position, 5 
it only contributes to the monopole 8oo(k) in Eq. (TToT l. and the 
monopole vanishes by definition at k = 0. Without loss of 
generality, we can express n g {z) as 



tg (z) = n g W(z) = 



V„ 



where the surveyed volume is 



W(z) , (19) 



V s = 4tt / dr r 2 W = 4tt / dz — W , (20) 

and JVtot is the total number of galaxies measured in the sur- 
vey. Equation f!9[ defines the survey window function W, 



5 By definition, the mean n g (z) is a quantity in a homogeneous universe, 
and the residual part S g (z, x) is a fluctuation around the mean. In obser- 
vation, the mean is independent of angle by construction. The two means 
should agree, provided that n 9 (x) is averaged over a large volume at a 
given redshift. 
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also known as radial selection function. It is related to the 
redshift distribution as 



VJz) 



H 



V r (r) 



4w r 



W(z) 



47r r A 
N^tH 



n g (z) , (21) 



where H(z) is the Hubble parameter and the normalization 
convention is 1 = J dz V z {z) — J dr r 2 V r {r). In principle, 
W(z) could be generalized to include an angular selection 
function. In what follows however, we will consider galaxy 
surveys with full-sky coverage and uniform angular selection 
function for simplicity. The galaxy number density at position 
x thus is 



n g (x) = h g W(z) [1 + S g (z,x)} , 
and its two-point correlation function reads 



(22) 



<n g (xi)n g (x 2 )> = n 2 ,W{z l )W{z 2 ) l + ^(x 2 -Xi) 



+n g W(z 1 )5 D (x 2 -x 1 ) 



(23) 



where £, g (x 2 — xi) is the Fourier transform of the noise-free 
galaxy power spectrum. We have assumed that the galaxies 
are an (inhomogeneous) Poisson sampling of 1 + S g (x) to de- 
rive the shot-noise term ll36ll . 

Before we proceed further, we introduce the transfer func- 
tions Tx (k, z) of random perturbation variables X that fur- 
ther simplify the spherical Fourier decomposition by separat- 
ing radial (time) dependence from angular dependence. For 
the galaxy fluctuation, we have 



5 g (k, z) = T g (k, z) ip v (k) + e(k, z) , 



(24) 



where the power spectrum of the comoving curvature 
^%S.k) — k 3 P V v(k) /2tt 2 at initial epoch is a nearly scale- 
invariant and the transfer function is independent of angle. 
The comoving curvature (f v (k) is often denoted as £(k) in 
literature. Arising from the discrete distribution, e(k, z) is 
a residual Poisson noise that is uncorrelated with S g (k, z), 
with power spectrum (e(k, z)e(k' , z)) — (27r) 3 <5 D (k — 
k')/(h g W(z)). 

In case that time evolution is related to radial coordinates, it 
is more natural to use Eq. ( [TBI than Eq. ( fT3l >, and so is it to use 
Eq. (TTTb than Eq. (TBI for computing spherical Fourier modes 
and their spherical power spectrum, respectively. On inserting 
Eq. (l22l into Eq. ( fT6] l and substituting the transfer function T g 
and survey selection function W, the spherical Fourier mode 
simplifies to 



Slm{k) 



dlnk'k 



2vr 2 



d 2 k! ^(k')y4(k') Mi(k' \k) 

(25) 



where ei m (k) is the spherical Fourier transform of the residual 
noise field h g W(z)e(x), and the spherical multipole function 
Mi(k' , k) is defined as 



dr r 2 W(r) ji{kr) j t (kr) T g (k,r) , 

(26) 



Mi(k,k) = kJ- 



where its dimension is [A4i(k, k)] — L 2 . After some simpli- 
fication, the spherical power spectrum in Eq. (TTTt eventually 
reads 

Si(k,k') = Ann 2 [dink A 2 (k)Mi(k, k)Mi{k, k') 



9 I — <p v 

2kk' 



TV 



+ n g / dr r 2 W{r)ji{kr)ji{k'r) .(27) 



o 



The second-term in the right-hand side is the shot-noise con- 
tribution. Using the Limber approximation (see Sec. IIII O . 
the shot-noise power spectrum can be rewritten as 



2kk' f°° 

M(k,k') = h g / drr 2 W{r)ji{kr)ji{k'r) 



fi g W{v/k)5 u {k-k') , 



(28) 



where v = I + 1/2. For a homogeneous and isotropic galaxy 
population with constant comoving number density n g = n g 
and power spectrum (S g (k)S*(k')) = (2n) 3 S D (k-k')P g (k), 
the spherical power spectrum Eq. (|27t yields the well-known 
relation 



5 i (fc) = n 2 9 P s (fc) + 



(29) 



The angular multipole I controls the transverse wavenumber 
k± but, since the amplitude of the wavevector k = \k\ is al- 
ready set in Si(k), the spherical power spectrum must be in- 
dependent of /. In practice, however, W(r) is always different 
than unity so that Si(k, k) is never independent of the multi- 
pole /. Note also that, with the galaxy number density defined 
as Eq. d22"i >. both the spherical power spectra Si and Mi have 
dimensions of L~ 2 . We will henceforth assume that V s is ac- 
curately known and work with the normalized galaxy number 
density n g (x)/n g , instead of n s (x). In this case, the spheri- 
cal power spectrum is given by the right-hand side of Eq. (l27t 
divided by h 2 . 

Before we close this section and apply the spherical Fourier 
analysis to the general relativistic description of galaxy clus- 
tering, we discuss our assumption for survey geometry and 
other approaches to handling the survey window function. 
Our spherical Fourier decomposition assumes that the full sky 
is available for measuring the galaxy number density field 
n g (x). For surveys with an incomplete sky coverage, an- 
gular multipoles with characteristic scale larger than the sky 
coverage are simply unavailable, while angular multipoles on 
smaller scales remain unaffected. 

In literature, the radial boundary condition is often imposed 
by choosing discrete wavenumbers k ni , in which the range of 
integer ni depends on angular multipole I (e.g., 1261 l27ll3 ill ). 
Since the galaxy number density can be measured only within 
the survey region, the number density field outside the sur- 
vey volume is set zero in those approaches. Hence, this sit- 
uation is equivalent to solving the Poisson equation within 
the survey area and the Laplace equation outside the sur- 
vey ll26ll . A unique solution can be singled out on imposing 
Dirichlet or Neumann conditions at the boundary (or a lin- 
ear combination thereof) |27|]. Therefore, despite the same 
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number density measured within the survey region, the result- 
ing potential and velocity fields are different due to the non- 
locality of their relation to the density field, depending on the 
choice of boundary conditions. As we are interested in de- 
scribing the observed galaxy number density n° bs (x), rather 
than the other derived quantities such as velocity or potential, 
no boundary condition need be imposed, and the results for 
n° bs (x) are identical in all approaches, whenever our contin- 
uous wavenumbers equal discrete ones or the survey volume 
becomes infinite. 



physical volume can be obtained by tracing the photon path 
backward (see iH [HI for details in deriving dV p hy and com- 
puting the volume effect). 

Another source of fluctuations in galaxy clustering is the 
source effect describing the distortion associated with the 
physical galaxy number density n v g hy . The mean galaxy num- 
ber density is obtained by averaging the observed galaxy num- 
ber density at the observed redshift, which differs from the 
proper time at the galaxy rest frame. The contribution of this 
source effect is proportional to the evolution bias factor of the 
galaxy number density 



III. GALAXY CLUSTERING IN GENERAL RELATIVITY 

Here we present the full relativistic description of galaxy 
clustering in Sec. IIII Al Weight functions are derived to facil- 
itate the full calculations of the spherical galaxy power spec- 
trum in Sec. IIII B I and simplified Limber formulas are pre- 
sented in Sec. IIII Cl to provide physical insight. 



A. General Relativistic Description of Galaxy Clustering 

Since the standard approach to modeling galaxy cluster- 
ing is based on the Newtonian framework, it naturally breaks 
down on large scales, where relativistic effect becomes sig- 
nificant. Recently, the fully relativistic description of galaxy 
clustering was formulated resolving gauge issues be- 

tween observable and unobservable quantities. The key find- 
ing in the theoretical development is that unobservable quanti- 
ties such as the true redshift and the unlensed galaxy position 
are gauge-dependent. Therefore, we need to construct galaxy 
clustering statistics based on observables quantities like the 
galaxy redshift and the position on the sky. In doing so, ef- 
fects such as redshift-space distortions and lensing magnifica- 
tion are naturally incorporated. Hence, this approach provides 
a unified description of various effects in galaxy clustering 

The relativistic description of galaxy clustering can be de- 
rived from the fact that the number dN° hs of observed galax- 
ies in a small volume is conserved: 



din? 



dN° hs 



n° bs dV ohs = nfvdV, 



phy 



(30) 



and this equation defines the observed galaxy number density 



-.obs 



with the observed volume element 



dV obs (z) = 



r 2 (z) 
H{z) 



sin 9 dz d6 . 



(31) 



described by the observed galaxy position (8, <fi) on the sky 
and the observed redshift z. The "observed" volume element 
dVcbs is different from the "physical" volume dV p hy occu- 
pied by the observed galaxies on the sky, and the distortion 
between these two volume elements, so called the volume ef- 
fect, gives rise to contributions to the observed galaxy fluctu- 
ation. While the dominant contribution in the volume effect is 
the redshift-space distortion and the gravitational lensing ef- 
fect, there exist other subtle relativistic contributions, and the 



dln(l + z) ' 



(32) 



where the factor three in Eq. 



appears because we use 



the comoving number density n g . Other contributions of 
the source effect can arise, depending on how we define the 
galaxy sample at hand (see ll9l- [Tll[T7ll ). In typical galaxy sur- 
veys, galaxy samples are defined with the observed luminos- 
ity inferred from the observed flux, which is different from the 
intrinsic luminosity. So the source effect associated with it is 
then proportional to the luminosity function slope 



P 



-0.4 



d log n g 
d log L 



(33) 



where the luminosity L is computed by using the observed 
flux and the observed redshift. 

Putting all these effects together and accounting for the rel- 
ativistic contributions, the general relativistic description of 
the observed galaxy fluctuation can be written, to the linear 
order in perturbation, as Jgl Holl 



jobs 



b S v — e Sz v — 5p 6T>l + a x + 2 ip x + V 



where the linear bias factor is b, the comoving gauge mat- 
ter density is S v , the lapse term Sz in the observed redshift is 
defined as 1 + z = (1 + z)(l + Sz), the dimensionless fluc- 
tuation in the luminosity distance is 6T>l, the gauge-invariant 
temporal and spatial metric perturbations are a x and <p x , the 
line-of-sight velocity is V, and the gauge-invariant radial dis- 
placement and lensing convergence are S1Z and JC. The sub- 
scripts x an d v represent that the gauge-invariant quantities 
with the corresponding subscript are identical to those quan- 
tities evaluated in the conformal Newtonian gauge (\ = 0) 
or in the dark matter comoving gauge (v — 0), in which the 
shear seen by the normal observer vanishes (\ — 0) or the 
off-diagonal component of the energy momentum tensor van- 
ishes (v = 0), respectively. Note that a x an d ip x correspond 
to the Bardeen's variables <&a and <&h 113711 and we have as- 
sumed no vector and tensor contributions in Eq. ((34). We refer 
the reader to JH [Tol [l7| for details in the derivation (see also 

MM)- 

Equation d34l is written in a physically transparent way. 
The observed galaxy fluctuation is modulated by the matter 
density S v , and there exist additional contributions from the 
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source and the volume effects. The source effect is composed 
of e 5z v and 5p 5T>l, and the other terms in Eq. ( l34l come 
from the volume effect. Of the volume effect, the last three 
terms in the first line of Eq. (l34l defines the rest frame of 
the observed galaxies, and the remaining terms in the second 
line describe the distortion of the comoving volume element 
in Eq. ( BIT ). It is noted that the term 3 5z x accounts for the 
distortion in the comoving volume factor a 3 and in Eq. (l3TT l 
the conversion between the comoving and the physical vol- 
ume elements is based on the observed redshift. 

In order to compute the observed galaxy fluctuation S° hs , 
we need to evaluate the individual gauge-invariant variables 
in Eq. d34l ). While Eq. ( [34-b is arranged in terms of gauge- 
invariant variables to explicitly ensure the gauge-invariance of 
Sg hs , it can be computed with any choice of gauge conditions, 
but preferably with the most convenient for computation. Us- 
ing the Einstein equations, we have the following relation for 
a flat ACDM universe JSSI, 



fx = 



3#o 2 



a. 



V — «y = 



2 ak 2 
1 A' U -f A 



(35) 
(36) 



where the comoving gauge matter density is obtained upon 
setting 6 m = S v . and the logarithmic growth rate is / = 
d In S m /d In a. Other gauge-invariant variables in Eq. (1341 can 
be readily expressed in terms of 5 m , v, and (f> as lfl7ll 



_ d f d 3 k —v 



dr J (2tt) 3 k 
5z x = V + cj>+ ( dr 24>' , 



5z v = Sz x + 



d 3 k Uv 



(2tt) 3 k 
611 = - 5 -^-J\f2cb, 



df [ r —^- | V 2 



511 



ST>l = h Sz x + (f> — JC , 

_ H ±(sl K \_ v _l_±z l 1 + zdV 

dz\n J n 

1 + zdH . 

—oz x H — — oz x 

x H dz 



H dr 



(37) 
(38) 
(39) 
(40) 
(41) 
(42) 

(43) 



We also note that the total derivative in Eq. d43l is along the 
past light cone: dr = dz/ H and dr — dr — dr. 



B. Fourier Decomposition and Weight Function 

Now that all the contributions in Eq. d34T > can be expressed 
in terms of the matter density S m and its transfer function 
T m (k,r), we define the weight functions M^(x, k) for each 
contribution in Eq. ( [34-b through the following decomposition 



of the observed galaxy fluctuation: 
d 3 k 



E 

i 

+ c(x) , 



dr 



(2tt) 3 



Wi(5t,k)T m {k,r) tp v (k)e 



ik-x 



(44) 



where the observed galaxy position is x = [r(z),x] and 
the line-of-sight position is x = (f, x) in spherical coor- 
dinates (x//x). The dimension of the weight functions is 
[Wi(x, k)] = and e(x) is the residual shot-noise field. 

Categorically, each contribution in Eq. (l34t can be classi- 
fied in three different types according to its dependence on 
the angle x and the line-of-sight distance f. The first and sim- 
plest type includes all the contributions such as S v , a x , and 
(p x that are independent of x and f. For example, the matter 
density in Eq. d34l ) is 



<5„(x) = 5 m (x) 



d 3 k 

J2^f 



T m (k, r) (p v (k) e 



ik-x 



(45) 



and hence the weight function for the matter density takes the 
simplest form 



W(x,k) ^S D (r-f). 
Similarly, the weight function for the potential <p x 

W(i,k) = S D (r-r) 3 -^^. 



(46) 



(47) 



Fluctuations of the second type depend on the observed an- 
gle, yet are independent of the line-of-sight distance. This 
applies to the line-of-sight velocity contribution V (and func- 
tions thereof such as Sz x , STZ, 5T>l), 



tti \ d f rf 3 k f S m 



(48) 



of which the weight function is 



W(x,k) = S D (i 



r '\k kdr 



(49) 



The derivative operator acts on the radial component of the 
plane wave in Eq. ( f44b . 

Finally, the third type covers the contributions along the 
line-of-sight direction that are independent of the observed 
angle. This category includes numerous line-of-sight inte- 
grals in Eqs. d3~8l~d42l. including the Sachs-Wolfe effect and 
the weak lensing convergence. For example, the Sachs-Wolfe 
contribution in Sz x in Eq. ( 1381 can be decomposed as 



r 

/ df (p 
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d 3 k 

(2tt) 3 



= / df 
Jo 



d 3 k 

(2tt) 3 



3^ H(f-l) 
2 m fc^ 



(50) 

T m (k,f) <p v (k) e lk -*, 
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TABLE I: Weight functions of the observed galaxy fluctuation in Eq. l !44t . All the perturbation variables in Eq. l !34t are scaled with the matter 
density in the dark-matter comoving gauge by using the Einstein equations, and its relations to the matter density are defined as the weight 
functions that can be used for computing the spherical multipole function in Eq. l |54t and the spherical power spectrum in Eq. l |57| l. Time 
dependent quantities a, H, and / are evaluated at the observed redshift z, and those with tilde depend on the line-of-sight distance f (or the 
corresponding redshift z < z). 

gauge-invariant quantity weight function Wi (r, r, k) 

Jo S D (r-f) 

-* D (r-f)3fa* 

**- * D <r - f) [(¥£0 + ^ - ^] + 3Hfa m M^ 

S D (r-f) \, H! 8 \ , 3i?gn m l 3ffg n H(f-l) 3Jfg Sl m 
r Hr [\ k kdr) 2 ~ak? \ Ur " m P r iF 

~H£ S D (r -)[(2^f)(M4) + El^ (/ ^ f ) + /^] _ SHSn m (l - M) 

/c 5|in mK/ + i)(z^f) i 

«>* ^ r) [(1 - (^) + ^ (2 - £)] + ^ (l _ _ £ _ Mi+I> ( ^)] 



and, therefore, the resulting weight function is 



(51) 



and the spherical multipole function of the observed galaxy 
fluctuation is simply 



(55) 



where H and / depend on f and we used <fi' — H(f>(f — 1). 
Another example of this type is the weak lensing convergence 
in Eq. gTJl 



K 



(If 



rr 
ri 3 k 
J2nf 



d 3 k 

r-f\ 3H$ 



<£(k, f) e ik i 



25fc 2 



xT m (fc,f) ^(k) VV 1 "* 



(52) 



The weight functions depend on position r, wavevector k, and 
line-of-sight distance f, and its angular dependence can be 
removed by using the plane wave expansion, i.e., 



W l (r,f,k) = ^n m l(l + l)f r -^j^, (53) 



where the expansion factor a depends on f and the angular 
part of the plane wave satisfies V 2 Y/ m (x) = — l(l+l)Yi m (x). 
For convenience, Table [Q summarizes all the weight functions. 

With the full Fourier decomposition of the observed galaxy 
fluctuation, the spherical multipole function in Eq. d26l > must 
be generalized as follows to take into proper account each 
component in Eq. d44l . 



The spherical power spectrum of the observed (normalized) 
galaxy fluctuation is then 

Si(k, k') = Si(k, k') +Ni{k, k') , (56) 

where Si(k, k') is the noise-free galaxy power spectrum, 

Si{k,k')=4n f dink Al v (k) Mi(k,k) Mi(k,k') , (57) 



and the noise component is 

Okh' roc 



's Jo 



drr 2 W(r)ji(kr)ji(k'r), (58) 



where Si is identical to Eq. d27b except for a factor of n 2 . Its 
dimension thus is L 4 . 

The weight functions defined here differ from those in H, 
as we scale all the perturbation variables with the matter den- 
sity in the comoving gauge. This approach is better suited 
to numerical computation, but its applicability is limited to a 
flat universe with a pressureless component, including a cos- 
mological constant and baryons on large scales. For different 
cosmological models or modified gravity theories, the relation 
among S v ,a x ,ip x , and v x is in general different [1421 - 44411 from 
what we assumed for a pressureless medium, and one should 
use the weight functions in [9] with the corresponding transfer 
functions. 



M\(k,k) 



10 roo 



dr r 2 W{r) ji(kr) (54) 



x / drWl{r,r,k)T m {k,r) j,(fcr) , 
Jo 



C. Limber Approximation 

Before we present numerical calculations of the spherical 
power spectrum in Sec. |IV] we adopt the Limber approxi- 



8 



mation to derive analytic formulas for the spherical multi- 
pole functions and the spherical power spectra. The Limber 
approximation (J45l-S3l) relies on the fact that the spherical 
Bessel function ji (x) peaks ata;«^ = Z + l/2 and oscillates 
rapidly for x > v. For large values of angular multipole I, in- 
tegrals of product of a smooth function f(x) times ji(x) can 
be approximated by 



(59) 



dx f(x) ji{x) ~ J — f(v) [1 + 0{l/v 2 )} 



With aid of the orthogonality relation for the spherical Bessel 
functions, this result can also be written as [481 



dr r 2 f(r)ji(kr)ji(k'r) 



S D (k-k') 
k 2 



(60) 



xf(u/k) [l + 0(l/v 2 )] 



In the case of the matter density, for instance, we obtain the 
spherical multipole function 




Mf(k,k) 



w W{v/k)T m 



k - 

k 'k 



8 D (k-k), (61) 



and the spherical power spectrum 



Sf(k,k') « W(v/k) 2 P„ 



k - 

k 'k 



S D (k-k') 



(62) 



Similarly, on using the Limber approximation, the spherical 
power spectrum Sf(k,k') of the potential perturbation can 
be readily computed owing to the simple dependence of the 
spherical multipole function on the weight function. The noise 
contribution to the spherical power spectrum falls in the same 
category, and the Limber approximation gives 

Ni(k,k')^ —W{v/k)5 D {k-k') . (63) 
n g 

In addition, the Limber approximation can greatly simplify 
the computation of the projected quantities such as the Sachs- 
Wolfe contribution and the gravitational lensing contribution. 
These quantities contain two line-of-sight integrations in the 
spherical multipole function in Eq. (l54l . each of which in- 
volves the spherical Bessel function. With the Limber approx- 
imation, the spherical multipole function of the gravitational 
lensing is 



Mf{k,k) 



w{v/k) n rn r m {k, v /k) 



xl(l + l) 



k-k 
dk 2 k 2 



(64) 



FIG. 1: Survey window functions VV and redshift distributions P z . 
Upper panel: Gaussian window function (solid) in Eq. d65t and top- 
hat window function (dotted) in Eq. J66l l. Bottom panel: Correspond- 
ing redshift distributions obtained by using Eq. J2U . The character- 
istic scales ro of the surveys are set equal to the comoving distance 
to 2 = 1. The comoving number density is related to the survey 
window function as in Eq. dl9t . 



The spherical multipole function of the line-of-sight veloc- 
ity or the redshift-space distortion involves the derivatives of 
the spherical Bessel function, since their weight function acts 
on the radial component as in Eq. (HSl l. Since derivatives of 
the spherical Bessel function ji(x) are linear combinations 
of spherical Bessel functions at different angular multipoles, 
they do not form an orthogonal basis. In these cases, we find 
that the Limber approximation becomes less accurate, and 
next-leading corrections are required to enhance the accuracy. 
Despite this shortcoming in a few cases, the Limber approx- 
imation provides nonetheless a simple and physical explana- 
tion of the complicated spherical power spectrum. Therefore, 
we will frequently rely on it when we discuss the physical in- 
terpretation of the full numerical results in Sec. [IV] 



IV. OBSERVED SPHERICAL POWER SPECTRUM 



where a is the expansion factor at f = v/k. The spheri- 
cal power spectrum Sf"(k, k') of the gravitational lensing is 
then obtained by integrating the spherical multipole function 
over k, instead of a quintuple integration over k and two pairs 
(r,f). 



Here we present the observed spherical power spectrum in 
Sec. IIV Al and compare to the corresponding flat-sky three- 
dimensional power spectrum. The measurement uncertainties 
associated with the spherical power spectrum are presented in 
Sec. UVBl 
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wavenumber k (hMpc" 1 ) wavenumber k (hMpc -1 ) wavenumber k (hMpc -1 ) 



FIG. 2: Spherical power spectra with Gaussian window function. In each panel, various curves with labels show the spherical power spectra 
of representative components of the full general relativistic formula S° hs in Eq. i34i , which are further decomposed in Eqs. (137) — ll43l. The 
time evolution of the transfer and the weight functions is ignored, and their values are set equal to those at z — 0. Each column plots the same 
spherical power spectra at angular multipoles, ranging from I — 2 (top rows) to / = 10 (bottom rows), illustrating their angular dependence. 
Vertical lines mark the characteristic wavenumber k c = l/ro, below which spherical power spectra are suppressed due to the constraint k > 
k± ~ k c . Gray curves show the three-dimensional power spectra at z = 0, if there exist well-defined correspondences. For example, spherical 
power spectra for the density (solid) and the redshift-space distortion (dashed) in the first column reduce to the corresponding three-dimensional 
power spectra on small scales. The noise power spectrum (dotted) is shown in the first column, assuming n g = 1(T 4 (h^Mpc)- 3 . Compared 
to the contributions in the first two columns, the projected quantities in the third column are negligible, except the gravitational lensing 
convergence. The turn-over of the spherical power spectra on large scales reflects the limit set by survey depth ~ a few times l/ro. These 
large-scale modes can be recovered by deconvolving the survey window function. 



A. Spherical Power Spectrum 



For definiteness we consider all-sky galaxy surveys with 
two different survey window functions. Figure Q] plots two 
survey window functions W(r) and their corresponding red- 
shift distribution V z {z). We consider two specific window 
functions that approximately represent a flux-limited survey 
(solid) and a volume-limited survey (dotted), if the comoving 
number density n g (z) is constant. Survey window functions 
(top panel) are related to galaxy redshift distributions (bottom 
panel) via Eqs. ( fT9b and d2lV and hence to the evolution bias 
factor in Eq. (I32t . The former window function is modeled as 



a Gaussian (solid) 

W(r) = exp 



(65) 



and the latter window function is modeled as a top-hat (dotted) 

W(r) = {;> r ^°< (66) 

where the characteristic distance or the survey depth is ro = 
2354/i -1 Mpc set equal to the comoving distance to z = 1. In 
both cases, the redshift distribution rises with redshift as more 
volume becomes available, but it drops sharply at ro for the 
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volume-limited survey, while the tail of the flux-limited sur- 
vey extends to higher redshift z > 1 . By defining the spherical 
Fourier modes through Eq. ([TBI , we are implicitly assuming 
that the galaxy overdensity <5 ff (x) is pixelized. Therefore, the 
survey window function W(r) should in principle include a 
pixel window function. For simplicity, however, we will ig- 
nore this complication and assume that the data has been de- 
convolved. 

Figure [2] illustrates the spherical power spectra of repre- 
sentative components of the full general relativistic formula 
gobs j n gq d34] ^ sxich as the matter density S m , the line-of- 
sight velocity V, the gravitational potential (j>, and so on (see 
Eqs. l37l — l43l ). The Gaussian window function is adopted 
in computing the spherical power spectra. For simplicity, 
we first assume that the matter transfer function T m (k, z) is 
independent of redshift (i.e., of the line-of-sight distance f) 
and the redshift-dependent quantities in the weight functions 
Wl(r, f,k) such as a, H, and / are set equal to those evalu- 
ated at z — when computing the spherical multipole func- 
tions A4](k, k) in Eq. ( 154-b . The spherical multipole function 
is then integrated over k to obtain the spherical power spec- 
trum Si (k) by using Eq. ( l57l i. The full spherical galaxy power 
spectrum Sf hs (k) is the sum of the auto and the cross power 
spectra of all the components in Eq. ( |34j >, which can be readily 
deduced from Fig.|2](see also BlOll ). 

The first column shows the spherical power spectra Si (k) of 
the matter density S m (solid) and the redshift-space distortion 
(dashed) as a function of wavenumber k and angular 
multipole I. Each row shows the same spherical power spec- 
tra but at different angular multipoles, ranging from I — 2 
(top) to I = 10 (bottom). Since the matter density power 
spectrum P m (k) is isotropic, the spherical power spectrum 
S[(k) (solid) of the matter density is independent of angu- 
lar multipole, and it reduces to the three-dimensional power 
spectrum P m (k) — Sf(k) (gray solid). Given angular multi- 
pole I in each row, the characteristic wavenumber k c = I /ro 
is shown as the vertical lines, below which the power is 
highly suppressed and well beyond which the equality of two 
power spectra holds. Technically, the power is suppressed at 
k <C k c due to the spherical Bessel function ji(kr) ~ at 
kr <C I, and the physical interpretation is that larger scale 
modes k <C k c ~ k± cannot be measured given the survey 
depth ro. The match between Sf(k) and P m {k) could be ex- 
tended to larger scales, if more survey volume is available than 
what is assumed here. 

For the Gaussian survey window function we adopted here, 
there exists an analytic formula for the spherical multipole 
function Aif(k',k). Using Eq. (6.633) in Gradshteyn and 
Ryzhik ll49tl and setting A = and fj, = v = I, a useful 
identity can be derived as 



dx x 2 e~ ax ' 



2 a 



(67) 



where the modified Bessel function is I v (z) = (—i) v J v {iz). 
If we again ignore the time evolution of the transfer function 
T m (k,z) ~ Tm(k), the spherical multipole function of the 



matter density is |2' 



M^k\k) = J^T m (k')^ 2 



and in the limit r^kk' S> 21(1 + 1) it reduces to 



jk-k') 2 



2 
(68) 



, (69) 



where the Gaussian wavepacket in the square bracket con- 
verges towards the Dirac distribution in the limit of infinite 
volume (ro — > oo). 

At sufficiently high k, where Aif (k ; , k) is well described 
by Eq. (l69l . the spherical power spectrum in Eq. (l57t becomes 



Sf(k,k) 



dkP m (k) 
ro 



ro 



-in 



Pm{k) 



2V2vr 



and we have 



Sf(k)=Sf(k,k) 



ro 



2V27T 



(70) 



(71) 



Equation d70b reduces to Eq. ( fT5l l, if the volume is infinite. 
Since numerical integrations of the spherical power spectrum 
in Eq. d57| > require fine sampling and prone to numerical er- 
rors due to the highly oscillating nature of the spherical Bessel 
functions in Mi(k, k), we tested our numerical calculations 
against the analytic solution, and both results agree to high 
precision. 

Next, the dashed curves in the first column represent the 
redshift-space distortion contribution ^j^^r- Since there is 
only one line-of-sight direction in the plane parallel limit, 
the redshift-distortion contribution is conveniently expressed 
in terms of cosine angle fi^ = x • k between the line-of- 
sight x and the wavevector k, and its three-dimensional power 
spectrum is P z (k, /ifc) = f 2 fi^,P m (k). Thin dashed curves 
show the spherical power spectrum Sf(k) of the redshift- 
space distortion, while the gray dashed curves in each row 
show P z (k,Hk = 1) = f 2 Pm(k) for comparison. As the an- 
gular multipole becomes higher at a given k = |k|, more an- 
gular modes are emphasized (or equivalently lower /i^), and 
the spherical power spectrum Sf(k) decreases. At k ^> k c , 
the spherical power spectrum reduces to the three-dimensional 
power spectrum P z (k,fj,k) with its angular dependence en- 
coded in angular multipole in each row. 

The second column plots the spherical power spectra of the 
dominant relativistic corrections to the observed galaxy fluc- 
tuation S° hs , such as the line-of-sight velocity V (solid), the 
gravitational potential <f> (dashed), and their linear combina- 
tions in Eqs. d37|i— d43l. Since the velocity and the potential 
contributions scale as (k/H) and (k/Ti) 2 relative to the mat- 
ter density 6 m , these relativistic contributions are comparable 
only at the horizon scale k ~ H and constitute small cor- 
rections in most galaxy surveys (see IU7I1 for measuring these 
effects). 
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FIG. 3: Spherical power spectra with top-hat window function, in the same format as in Fig. [2] The general trend in the spherical power spectra 
are similar, but oscillations are imprinted in the spherical power spectra, due to the incomplete cancellation of the spherical Bessel function, 
imposed by the top-hat window function. The gray curves are identical to those in Fig. [2] 



The spherical power spectrum Sf(k) (dashed) of the gravi- 
tational potential <f> follows the three dimensional power spec- 
trum P<p{k) (gray dashed) at each row, whose amplitude 
quickly decreases at high k as P$(k) oc k~ 3 , but is in- 
dependent of angular multipole I. At a given small k, the 
amplitude Sf(k) becomes lower in the lower panels, sim- 
ply because k < k c (the vertical line in each row again il- 
lustrates the characteristic minimum wavenumber k c given 
angular multipole I). The angular dependence of the line- 
of-sight velocity V means that its spherical power spectrum 
SY (k) (solid) also depends on angular multipole I and it de- 
creases with increasing angular multipole (or smaller /x^) on 
all scales, in addition to the effect at k < k c . Since the 
three dimensional power spectrum of the line-of-sight veloc- 
ity is P v {k^ k ) = fi 2 k P v (k) = {Hf/k) 2 fi 2 k P m {k), we plot 
Py(fc,/ifc = 1) (gray solid) at each row for comparison, 
and again the spherical power spectrum S t ■ (fc) reduces to the 
three-dimensional power spectrum Pv(k, /ifc). Last, the dot- 
ted curves show the noise power spectrum Afi(k) for a galaxy 
sample with n g = 10~ 4 (/i _1 Mpc) -3 . For the Gaussian ra- 
dial selection function, the shot-noise contribution to Si (k, k') 



can be computed exactly by using Eq. d67| i as 

M(fc,fcO = ^VIFe-^( fc2 + fc ' 2 )/ i+| ^) , (72) 

and it asymptotically matches the Poisson power spectrum 
l/fig on small scales. 

The remaining two curves show the variants of the line- 
of-sight velocity and the gravitational potential, V/Hr (dot 
dashed) and (j)/T-Lr (dotted). They both arise due to the ra- 
dial distortion 811/ r. With the prefactor 1 /Wr in those terms, 
more weight is given to shorter distance 7-, and hence they 
peak at a slightly higher wavenumber k than their counter- 
parts without the prefactor, while the shift in the peak position 
depends on the adopted window function. 

The last column shows the spherical power spectra of the 
projected quantities along the line-of-sight direction, such as 
the gravitational lensing JC (solid), the integrated Sachs-Wolfe 
effect J Q r df (j)' (dashed), the Sachs-Wolfe effect fl dr <p/r 
(dot-dashed), and their linear combinations in Eqs. (137t— d43ll. 
Since these contributions are intrinsically angular and their 
variations are largely limited to the transverse direction, the 
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angular power spectrum C; is often used to characterize them, 
and there is no exact analogy to the usual three-dimensional 
pow er spectrum for the projected quantities (see, however, 
IdOIO . Apart from the gravitational lensing convergence, the 
spherical power spectra of the projected quantities are at least 
an order-of-magnitude smaller than the gravitational poten- 
tial <j> in the second column. Those projected quantities are 
isotropic and independent of angular multipole, with the sole 
exception of the gravitational lensing JC, whose additional an- 
gular dependence 1(1 + 1) arises due to the angular Lapla- 
cian operator. The spherical power spectrum of the gravita- 
tional lensing JC exactly corresponds to the 3D weak lensing 
lHHIHlHj]], where the same spherical Fourier analysis is ap- 
plied to background source galaxies to map the foreground 
matter distribution. While two more shear fields are available 
in weak lensing measurements, only the convergence field JC 
contributes to galaxy clustering. 

In Fig. [5J we show the spherical power spectra with the 
top-hat window function in Eq. ( f66b and discuss its effect 
on the spherical power spectra, compared to those in Fig. |2 
We have again ignored the time evolution of the transfer and 
the weight functions for simplicity. Three-dimensional power 
spectra shown as gray curves in Fig. |3]are identical to those in 
Fig. |2] The normalization of the spherical power spectra can 
be obtained in a similar manner by considering the spherical 
power spectrum of the matter density, 



S?(M) = / dkP m {k) 



2 kk 



dr r 2 ji(kr)ji(kr) 



and taking the limit k — > oo (see AppendixlAb 



lim Sf(k,k) = P m (k)^J- 

fc->oo TT \ IT 



(73) 



(74) 



Therefore, we plot the spherical power spectra in Fig. |3]as 

/ , \ -l 



S l (k)=S l (k,k)l^J^ 



(75) 



The spherical power spectra of the matter density (solid) 
and the redshift-space distortion (dashed) in the first column 
are nearly identical to their counterparts with the Gaussian 
window function. The prominent difference is the oscilla- 
tions in Si(k) at k < 27r/(ro/3), where most contributions 
to Si(k) at those low k arise near r ~ ro and the highly os- 
cillating spherical Bessel function is abruptly truncated at the 
survey boundary r = ro, causing incomplete cancellation of 
the oscillations. The noise power spectrum (dotted) is nearly 
constant under the top-hat survey window function. With full 
weight given W(r) = 1 near r < ro, however, the recovery 
of the three-dimensional power spectra is somewhat extended 
to lower k than in Fig. [2] 

A similar pattern is observed for the line-of-sight velocity 
(k) and the potential Sf(k) spherical power spectra in the 
second column. Due to additional suppression factor 1/fc 2 in 
the weight function for the potential <f>, more weight is given 



to large r at a given k for <f> than for V, which results in oscil- 
lations in Sf(k) even at k > 27r/(ro/3). The remaining two 
curves (dot-dashed and dotted) in the radial distortion follow 
the same trend but, compared to Fig. [2] their peak is some- 
what shifted to larger scale k due to the change in the window 
function. 

The spherical power spectra of the projected quantities in 
the third column show the reversed trend, smooth power on 
large scales and oscillations on small scales, compared to 
the spherical power spectra in the first column. The pro- 
jected quantities have additional integral in the spherical mul- 
tipole function A4i(k, k), and this line-of-sight integration 
cancels the oscillating part of the second spherical Bessel 
function, while the overall oscillation is determined by the 
first spherical Bessel function. The high-fc oscillation arises 
near r l/k, at which there is no suppression in the top-hat 
window function. Furthermore, the peak position is shifted to 
smaller scales, since the mean distance to the source galaxies 
is smaller with the top-hat window function. 

Having understood the key features in the spherical power 
spectra with the simplifying assumptions, we are now in a 
position to consider the full generality: time dependence of 
the transfer function and the weight function on the line-of- 
sight distance. In Fig.|4j we present the spherical power spec- 
tra with the Gaussian window function, explicitly accounting 
for their time dependence in the spherical multipole function. 
Overall, the spherical power spectra in Fig. |4]resemble those 
in Fig. |2] where the time dependence is neglected. The key 
difference in this case lies in the amplitude of the spherical 
power spectra, as the matter transfer function T m (k,r) de- 
creases in amplitude at r > (i.e., z > 0). 

The first column shows the spherical power spectra of 
the matter density (solid) and the redshift-space distortion 
(dashed). Compared to the fiat-sky matter power spectrum 
(gray solid) at z = 0, the spherical power spectrum Sf(k) 
is lower in amplitude, as it is obtained by averaging over a 
range of redshift in Fig.Q] For the same reason, the logarith- 
mic growth rate / is higher at z > 0, and hence the spherical 
power spectrum of the redshift-space distortion is higher in 
amplitude. The noise power spectrum (dotted) remains un- 
changed in Fig. |U since the time dependence of the mean 
number density is already accounted by the survey window 
function in Fig. [2 

In the second column, the spherical power spectra of the 
line-of-sight velocity, the gravitational potential, and their 
variants show little difference, compared to those in Fig. |2 
A small change in amplitude arises due to the shift in the peak 
line-of-sight distance. The velocity is nearly constant in red- 
shift, while the gravitational potential is somewhat larger at 
z > 0, as dark energy domination leads to the decay in the 
gravitational potential. This trend is further reflected in the 
third column, where the projected quantities are displayed. 
The spherical power spectra of the gravitational lensing (solid) 
and the gravitational potential (dot-dashed) are slightly larger 
in amplitude, while the integrated Sachs-Wolfe contributions 
(dashed and dotted) are smaller, since the potential is constant 
in the Einstein-de Sitter phase at z > 0. 
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FIG. 4: Spherical power spectra with Gaussian window function. The time evolution of the transfer and the weight functions over the line-of- 
sight distance is explicitly computed. Compared to the spherical power spectra without evolution in Fig. [2] the spherical power spectra remain 
largely unaffected in shape, while the amplitude is naturally lower than the three-dimensional power spectra (gray) at z — 0. Various curves 
are in the same format as in Fig. [2] 



B. Statistical uncertainties 

The covariance matrix of the spherical power spectrum 
Si(k,k) generally involves the power spectrum, bispectrum 
and trispectrum of the spherical Fourier modes Si m (k). Non- 
Gaussianity induced by nonlinear gravitational clustering or, 
possibly, already present in the initial conditions can gener- 
ate significant covariance through the bispectrum and trispec- 
trum. For simplicity, however, we will henceforth ignore these 
contributions and assume that 8i m (k) is simply Gaussian. 

We begin with the calculation of the covariance of the 
spherical power spectrum assuming an all-sky experiment. At 
fixed multipole I, Si (k, k) can be estimated by averaging over 



all to and a thin wavenumber interval Ak and by subtract- 
ing the shot-noise component. However, since the shot-noise 
contribution A/j is difficult to compute without precise knowl- 
edge of the radial selection function, we simply present the 
observed spherical power spectrum, i.e., without shot-noise 
subtraction. We thus define 

m Ak 

where the integration runs over the domain [k — Ak/2, k + 
Ak/2. The covariance of this (unbiased) band-power estima- 
tor is 



COT[£,(M),4(fc',*')] 



(21 + 1) (21' 



, ( Afc ) 2 



du / du' 



(|^W| 2 |W(w')| 2 >-(|^mW| 2 >(|WK)| 2 ) 

(77) 
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The assumption of Gaussian random fields considerably simplifies the evaluation of the trispectrum. A straightforward applica- 
tion of Wick's theorem yields 



Cav[Si(k,k),Sy(^,kf)] 



1 



(21 + if ^, (Afc) 

V 1 m.m' 



E 



i 



du \ du 

Afc J Afc 



Sl(u,u') +Afl(u,u') (5 m , m > +S m .-rn')Sli 



21 + 1 {Aky 



du / du' 

Afc " Afc 



St(u,u') +M(u,u') 



(78) 



In the flat-sky approximation, estimates from different band- 
powers are uncorrelated in the Gaussian limit, so that the co- 
variance of power spectrum estimators Cov[P(fcj), P(kj)\ is 
proportional to 2P(ki)P(kj)5k i ,k j - In the spherical Fourier 
decomposition, the Limber approximation shows that, for the 
matter (and potential) perturbations, band-power estimate of 
the spherical power spectrum are also uncorrelated. How- 
ever, contributions from, e.g., the line-of-sight velocity or z- 
distortions will introduce covariance between measurements 
from different band-powers. 

For shorthand convenience, let us now define 



Siik) = S t (k) + M{k) 



(79) 



If we restrict ourselves to a galaxy sample of constant comov- 



ing density n g = ft 
error on Si(k) scales as 



[A5,(fc)] 2 = 



N s /V s and Poisson white noise, the 



2/ + 1 



Si(k) + — 
n„ 



(80) 
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for infinitely narrow k bins. This expression is similar to 
the FKP expression [36], except that the number N(k)/2 of 
independent modes in a fc-shell is replaced by the number 
(21 + l)/2 of independent \klm) modes at fixed k and I, 
In general, however, the survey window function is not uni- 
form and I /fig must be replaced by the line-of-sight integral 
Afi (k). The cumulative signal-to-noise for a measurement of 
Si is then 



( Mk) 



FIG. 5: Signal-to-noise for a measurement of the spherical power 
spectrum of galaxies with constant bias b = 2 (see text). Results 
are shown as a function of the maximum multipole Z ma x- The curves 
show Eq. d85b for the tophat (TW) and Gaussian (GW) survey win- 
dow with several choices of ro. The filled symbols show S/N in the 
flat-sky approximation (Eq. 11851 ) for the Gaussian survey window 
with r = 2354 (circle), 1177 (square) and 589 h~ 1 Mpc (triangle). 
Note that their abscisse is arbitrary. 



(81) 



v l/3 



1=2 



21 + 1 



1/3 

where the factor of V s has been introduced so that the 
signal-to-noise has the appropriate units (i.e., dimensionless). 
The factor of 1/2 arises from the fact that i5(x) is real, i.e. 
3*m(k) — 8i- m (k). Note that the cumulative signal-to-noise 
for a measurement of, e.g., the spherical power induced by 
the line-of-sight velocity is given by the same formula with 
SY replacing Si in the numerator. 

In order to establish a correspondence with mode counting 
in the 3D flat case, let us consider the limit Si Afi. Let also 

1 /3 

L = V s and AL be the characteristic length and resolution 

of the survey. The maximum wavenumber and wavenumber This is equal (up to a factor of order unity) to the signal-to- 



*dk 
2^ 



Si(k) 



M(k)+Si(k) 



resolution thus are fc ma x = 2tt/AL and Afc = 2tt/L. The 
angular resolution of the survey is 8 ~ AL/L, which yields a 



maximum multipole l r 
signal-to-noise is 



1/6 ~ fc max /Afc. Therefore, the 



^max V max + 2) 



Ak 



(82) 



Afc 
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noise in the flat-sky limit, 

2 



x 2tt 



dk± fcj 



dkn 

~2T 



(83) 



Afc 2 



k — k ■ 

Afc 



/ "-max 

2 V Afc 



where the last equality assumes k m - m <C k max . When deal- 
ing with real galaxy data however, one must take into account 
the shot-noise and the survey window function. Ignoring con- 
tributions others than fluctuations in the matter density, the 
Limber approximation gives 



Kin 



yl/3 /max 

2 ^ 

2=2 



°°dk 
2^ 



21 + 1 



(84) 



n g W($)P m (k,%) 



In Fig. [5] this expression is evaluated as a function of Z max for 
the Gaussian (GW) and tophat (TW) survey window consid- 
ered in Sec. ^IV I (dotted and solid curves). The galaxy bias is 
assumed to be constant, 6 = 2, and the z-dependence of the 
matter power spectrum has been neglected, i.e., P m (k, r) = 
P m (k, 0). To exemplify the scaling (S/N) cx V}^ & , results 
are also shown for the Gaussian window with 1/8 and 1/64 the 
fiducial volume V s ps 55/i _1 Gpc. The signal-to-noise even- 
tually saturates at large values of / max because, for / ^S> 1, 
the Gaussian window W(v/k) ~ exp[— (l/kro) 2 ] only picks 
up (small-scale) measurements that are shot-noise dominated. 
For our choice of b, this occurs at k ~ 0.2/i _1 Mpc. As a re- 
sult, the signal-to-noise flattens out around (* max ~ fcro —400 
(r = 2354/i~ 1 Mpc) and -100 (r = 589/i _1 Mpc). The 
filled symbols represent the signal-to-noise in the flat-sky ap- 
proximation, 



— I = — I drr" 

N 7T 



dk k 2 



n g W(r)P m (k,r) 



10 Jo K l + h g W(r)P m (k,r), 

(85) 

computed for the Gaussian survey window with ro = 2354 
(circle), 1177 (square) and 589/i — 1 Mpc (triangle). The flat- 
sky estimates are consistent with the Limber-approximated 
(S/N)(< Z max ) m me limit l max 1. This provides sup- 
port for the validity of Eq. 1 



V. DISCUSSION 

We have performed an all-sky analysis of the galaxy power 
spectrum, accounting for all the relativistic effects in galaxy 
clustering. The spherical Fourier analysis has been well de- 
veloped in galaxy clustering jUE3l> while its application 
to galaxy clustering has been limited to the Kaiser formula 
115 211 . We have used the spherical Fourier analysis to analyze 
the full relativistic formula. The observed galaxy fluctuation 



is decomposed in terms of spherical harmonics and spheri- 
cal Bessel functions that are angular and radial eigenfunctions 
of the Helmholtz equation, providing a natural basis for the 
observer at origin to describe the galaxy clustering measure- 
ments on the observed sphere. 

In light of the recent development in the relativistic formu- 
lation of galaxy clustering ffl Holl . there exist numerous rela- 
tivistic effects in galaxy clustering, in addition to the standard 
redshift-space distortion effect. These relativistic effects be- 
come substantial on very large scales, and measurements of 
these large-scale modes inevitably invoke complications asso- 
ciated with the flat-sky approximation and the survey geome- 
tries. By using the spherical harmonics for its angular decom- 
position, the spherical power spectrum is independent of the 
validity of the flat-sky approximation, while it retains the ad- 
vantage of the standard Fourier analysis, namely, the simple 
and physically intuitive interpretation of the measurements in 
conjunction with the underlying matter distribution. 

We have computed the spherical Fourier power spectrum 
of the observed galaxy distribution, accounting for the rel- 
ativistic effects. Compared to the standard Newtonian de- 
scription, there exist additional contributions to the observed 
galaxy fluctuation, and these additional contributions can be 
categorized as the matter density fluctuation, the line-of-sight 
velocity contribution, the gravitational potential contribution, 
and the line-of-sight projection contribution IIL?T - [l6ll . The 
spherical power spectrum of the matter density is identical to 
the usual three-dimensional matter power spectrum Sf(k) ~ 
P m (k), regardless of its angular multipole I, as the matter 
power spectrum is isotropic. This correspondence greatly fa- 
cilitates the physical interpretation of the measurements. 

Since the line-of-sight velocity affects the observed dis- 
tance to the galaxies in redshift space, the velocity contri- 
bution in galaxy clustering is angular dependent. Therefore, 
the spherical power spectrum (fc) of the velocity contribu- 
tion is similar to its three-dimensional counterpart Py (fc, /Ufc), 
but its angular dependence is encoded as a function of angu- 
lar multipole I: At a given amplitude k of a wavevector k, 
higher angular multipoles represent larger transverse modes, 
or lower cosine angle fj,k = fcii/fc. The redshift-space dis- 
tortion effect arises from the spatial derivative of the line-of- 
sight velocity, and its spherical power spectrum Sf(k) follows 
the similar trend, as its three-dimesional power spectrum is 
P z {k^ k )=f 2 tiP m (k). 

The gravitational potential contribution to the observed 
galaxy fluctuation can be readily computed from the spher- 
ical Fourier decomposition. Similar to the case of the mat- 
ter density fluctuation, the gravitational potential power spec- 
trum P<j>(k) is isotropic, and the spherical power spectrum for 
the gravitational potential is identical to the three-dimensional 
power spectrum Sf{k) ~ P^{k), regardless of its angular 
multipole. Furthermore, while the gravitational potential con- 
tribution to the variance of galaxy clustering often diverges 
due to its scale-free nature, the spherical power spectrum is 
unaffected by this problem since it measures individual modes 
of fluctuations as in the traditional power spectrum analysis. 

The other contributions to the observed galaxy fluctuation, 
such as the gravitational lensing and the integrated Sachs- 
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Wolfe effects, arise from fluctuations along the line-of-sight 
direction. Since these projected quantities are intrinsically an- 
gular, it is difficult to handle their contribution in the standard 
power spectrum analysis. However, with its angular decom- 
position using spherical harmonics, the spherical Fourier anal- 
ysis can naturally implement the contribution of the projected 
quantities to galaxy clustering. With the sole exception of the 
gravitational lensing effect, we find that compared to the mat- 
ter density fluctuation, the contribution of the projected quan- 
tities are negligible, and this justifies the simplification of ig- 
noring the projected quantities in the power spectrum analysis 
ifiol [l7ll . Moreover, the spherical power spectrum Sf(k) of 
the gravitational lensing contribution in galaxy clustering is 
also known as 3D weak lensing 1I32I.I33I1 . where the informa- 
tion on radial distances to the background source galaxies is 
utilized to map the matter distribution in 3D, as opposed to the 
traditional 2D weak lensing. Our spherical Fourier analysis 
provides a complete and comprehensive description of galaxy 
clustering and its associated effects. 

We have derived the covariance matrix of the spherical 
power spectrum, assuming that the matter density fluctua- 
tion is the dominant contribution. The covariance matrix of 
the spherical power spectrum asymptotically matches that of 
the three-dimensional power spectrum on small scales. It is 
also shown that since the volume available for 

galaxy surveys at low redshift is relatively small, there are too 



few large-scale modes that are sensitive to the relativistic ef- 
fect in galaxy clustering. Therefore, it makes little difference 
in terms of measurement significance, if one chooses to em- 
bed the observed sphere in a cubic volume and performs the 
standard Fourier analysis, instead of performing the spherical 
Fourier analysis. However, the spherical Fourier analysis pre- 
sented in this paper provides a more natural way to analyze the 
full relativistic effects in galaxy clustering. Furthermore, it is 
shown [17] that the multi-tracer analysis [15311 with the shot- 
noise cancelling technique ll54ll can substantially enhance the 
measurement significance of the relativistic effects in galaxy 
clustering, in which case we expect that the spherical Fourier 
analysis becomes essential in describing the measurements on 
large scales. 
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and taking the limit (k — > oo). We first define the integrand 



F{k,k) 



2 kk 



7T Jo 

and then arrange Eq. (IAU as 



dr r 2 ji{kr)ji{kr) , (A2) 



/OO 
dk P m (k + k) [F(k, k + k)] 2 . (A3) 
-oo 



Since the integrand peaks around k ~ k, we take the limit 
(k — > oo) and expand the integrand F(k, k + k): 



lim Sf(k,k) P m (k) / dk 

k— >oo 
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where we performed a Gaussian integral. The asymptotic val- 
ues of the integrand are 
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and we obtain the normalization coefficient in Eq. d74l 



lim Sf(k,k) = P m (k) 
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Appendix A: Top-hat normalization 



Here we derive the normalization coefficient for the top-hat 
window function in Eq. (l74l . The normalization coefficient 
can be obtained by computing the spherical power spectrum 
for the matter density 



Sf(k,k) 





2kk r 
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Appendix B: Covariance of the spherical power spectrum 



Here we present details of the calculation of the covari- 
ance matrix of the spherical power spectrum estimator. We 
will omit the band-power averaging for the sake of concise- 
ness. Using Eq. (TToT i. the four-point correlator of the spherical 
Fourier modes reads 
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{S lm (k)5t m (k)5i, m ,(k')Sp m ,(k')) = 



|) k 2 k ,2 f[yd 3 x l ^ji(kr 1 )j l (kr 2 )j l (k'r 3 )j l {k'r 4 ) 
xy4(x 1 )y im (x 2 )y i f m ,(x 3 )^' ro '(x4)(5(x 1 ) ( 5(x 2 )5(x 3 )5(x4)) . (Bl) 



Assuming the galaxy overdensity field <5(x) follows Gaussian statistics, the four-point correlator in the right-hand side reduces to 



the sum of three products of two-point correlation functions, whose explicit expression is given by Eq. (|23l l divided by n g . The 
contribution that involves ((5 g (xi)£ g (x2))(<5 g (x3)(5 g (x 4 )) exactly cancels out the term (\8i m (k)\ 2 )(\8i' m i(k')\ 2 ) in the covariance 
matrix in Eq. ( TTTb . To compute the two other contributions, we must evaluate, e.g., 

J14 = -kk' fdr 1 r' 2 fdr 4 r\ ji(kri)jp (k'r 4 ) /d 2 xi fd 2 x 4 Y^x^Yy m > (x 4 ) 
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On expressing £ s (xi — x 4 ) as the Fourier transform 



x i W(n)W(r 4 ) [l + e ff (xi - x 4 )l + ^W(n)^(xi - x 4 ) 



,ik-(xi-x 4 



/d 3 k 
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and inserting the Rayleigh expansion in Eq. (0, the integrals over the angular variables simplify to 
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As a consequence, the contribution of £ 9 (xi — x 4 ) to J\ 4 becomes 
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Similarly, using <5 D (xi — x 4 ) = £ D (ri — r 4 )<5 (xi — x 4 )/r 2 and the orthonormality of the spherical harmonics, the contribution 
of the Poisson noise term to Ji 4 is 



dnr 2 /dr4raW(ri)j,(*ri)ii'(fc / r4)^(r 1 -r4) /d 3 iil£ l (ii)lW(*i) =M{k,k') S W S.. 



(B6) 



Ji 4 is the sum of these two contributions, i.e. J 44 = [Si(k, k') + Afi(k, k')]8w 8 mm i . Symmetry considerations show that this 
holds true for Jy, i ^ j. This leads to the desired result in Eq. (l78l >. 



